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Abstract.  The  information  preservation  (IP)  method  is  extended  in  this  paper  to  include  the  case  of  temperature 
variation.  The  IP  velocity  u  and  IP  energy  e  (or  IP  temperature  T)  are  introduced.  They  migrate  with  the 
simulated  molecules  and  are  governed  by  the  momentum  and  energy  conservation  equations.  The  method  is 
illustrated  on  the  example  of  rarefied  gas  between  parallel  plates  with  different  temperatures.  Comparison  of  the 
IP  simulation  result  is  made  with  experimental  data,  DSMC  calculation  and  solution  of  non-linear  Boltzmann 
equation. 


INTRODUCTION 

The  unavoidable  statistieal  scatter  that  is  present  in  most  DSMC  calculations  is  generally  regarded 
as  the  most  serious  practical  and  ideological  difficulty  with  the  method,  and  becomes  insurmountable 
for  low  Mach  number  flows  characterized  by  a  small  ratio  of  macroscopic  to  molecular  velocity.  For 
low-speed  flow,  useful  small  macroscopic  information  may  be  completely  lost  in  the  back  ground 
noises.  And  as  the  deviation  of  fluctuation  decreases  as  the  reciprocal  of  the  square  root  of  the  sample 
size,  the  signal  may  remain  submerged  even  for  rather  larger  sample  sizes  which  require  enormous 
computation  time.  In  order  to  solve  this  problem  and  make  use  of  the  collectivity  of  simulated 
molecules  in  the  DSMC  method,  the  information  preservation  (IP)  method  was  proposed.  The 
essence  of  the  technique  is  to  store  besides  the  ordinary  thermal  velocity  of  the  simulated  molecules, 
also  the  so  called  information  velocity,  that  is  to  record  the  collective  velocity  of  the  enormous  number 
of  real  molecules  that  the  simulated  molecule  represents.  Macroscopic  velocity  and  shear  stress  are 
obtained  from  the  time  or  ensemble  average  of  the  information  velocities  of  simulated  molecules  in  a 
cell  or  on  a  wall  element.  The  method  has  been  applied  to  imidirectional  flows  such  as  Couette  flow, 
Poiseuille  flow  and  Rayleigh  flow.  Excellent  agreement  is  obtained  with  exact  solutions  of 
continuum  and  free  molecular  flows  and  the  numerical  solutions  of  the  linearized  Boltzmann  equation 
in  transition  regime.  The  technique  is  being  developed  to  treat  the  two  dimensional  and  three 
dimensional  low  speed  rarefied  gas  flows,  i.e.  the  two  dimensional  flow  around  a  NACA  0012  air 
foil  and  the  three  dimensional  flow  in  a  micro-machined  membrane  particle  filter.  It  is  natural  to  extend 
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the  IP  technique,  enabling  it  to  deal  with  situation  in  which  temperature  varies  in  the  flow  field  not  due 
to  retardment  of  the  flow  but  due  to  the  transfer  of  heat  from  the  non-isothermal  boimdaries. 
Unfortunately  this  extension  is  not  a  trivial  one.  The  nontriviality  is  caused  by  the  fact  that  the  average 
energy  flux  of  a  monatomic  molecules  cross  a  surface  element  in  a  stationary  gas  is  IkTT^  (where 
is  the  particle  number  flux  =ncl  4,  here  n  is  the  number  density  of  the  flow  and  the  c  is  the  mean 
molecular  speed),  whereas  the  average  energy  of  a  monatomic  molecule  is  3kT  12,  so  that  a  simulated 
molecule  carrying  an  energy  of  amount  3kT/2  could  not  reexhibit  the  physically  real  energy  flux 
IkTF^.  Despite  this  difficulty  in  this  paper  an  attempt  is  made  to  enable  the  IP  method  to  include  the 
case  of  temperature  variation,  which  takes  place  in  some  flow  situation  of  micro-electro-mechanical- 
systems  (MEMS).  The  present  version  is  known  to  have  a  certain  error  owning  to  the  above  mentioned 
relationship  concerning  to  the  energy  amount  and  energy  flux.  The  procedure  is  outlined  in  the  next 
section  followed  by  the  illustration  on  the  example  of  temperature  and  density  distribution  of  rarefied 
gas  between  parallel  plates  with  different  temperatures.  The  results  are  given  for  rather  large 
temperature  difference  (ATIT^),  they  are  in  agreement  for  both  heat  transfer  and  density  and 
temperature  distribution  for  large  Kn  numbers.  As  for  small  Kn  numbers  IP  procedure  yields  density 
and  temperature  results  in  fairly  good  agreement  with  DSMC  simulation  but  gives  heat  transfer  results 
with  apparent  discrepancy  from  the  DSMC  simulation  results  for  small  Kn  numbers. 


OUTLINE  OF  THE  PROCEDURE  OF  IP  METHOD  FOR 
THE  CASE  OF  TEMPERATURE  VARIATION 


Any  IP  procedure  is  implemented  on  the  basis  of,  and  embedded  in  a  successfully  operated  DSMC 
code.  Except  the  usual  microscopic  (thermal)  velocity,  macroscopic  velocity,  density,  and  temperature 
introduced  in  the  DSMC  method,  the  IP  (information  preservation)  velocity  u  and  IP  energy  e  of  the 
simulated  molecule  are  introduced  as  well  as  IP  velocity  {/ ,  IP  density  p  and  IP  temperature  T  of  a  cell 
which  represent  the  macroscopic  properties  of  a  cell.  The  IP  velocity  u  and  IP  energy  e  migrate  with 
the  simulated  molecules  and  are  governed  by  the  momentum  and  energy  conservation  equations: 
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o-ldS 


ll|p  ^  = -IP  ■  + IP<^  ■  ■  >■ 


Ids 


(1) 


(2) 


where  conventional  notation  is  used  and  /  is  the  unit  normal  to  the  surface.  Meanwhile  the  IP  velocity 
and  temperature  of  a  cell  are  obtained  as  the  weighted  averages  of  the  IP  values  of  simulated  molecules 
in  the  cell.  The  macroscopic  density  of  a  cell  is  governed  by  the  mass  conservation  equation: 


fOP  =  pU-ldS 


(3) 


For  isothermal  flows  equation  (2)  can  be  omitted,  and  in  some  cases  a  can  be  simplified  with  normal 
pressure  alone  retained.  For  slow  non-isothermal  flows  equation  (2)  is  needed  and  the  Burnett  terms 
together  with  viscous  terms  in  o  need  to  be  incorporated.  The  concrete  procedure  and  formulation  will 
be  specified  for  each  example  problem. 
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RAREFIED  GAS  BETWEEN  PARALLEL  PLATES 
WITH  DIFFERENT  TEMPERATURES 

Rarefied  gas  is  confined  between  two  parallel  plates  at  rest  (one  aligned  on  y=-D/2,  the  other  on 
y=D/2)  with  different  temperatures  -  AT  and  1^  +  AT .  The  Knudsen  number  Kn  is  defined  as  D , 
where  Ag  is  the  mean  free  path  at  y=0.  The  study  of  heat  transfer,  density  and  temperature  distribution 
in  the  gas  for  a  variety  of  Knudsen  numbers  have  been  pursued  by  many  authors.  The  non-isothermal 
IP  procedure  is  illustrated  on  this  problem.  For  simplicity  the  monatomic  molecules  are  considered. 
The  movement  of  the  simulated  molecules,  the  encounter  of  them  with  body  surface  and  the  collision 
between  them  are  simulated  by  the  DSMC  method.  After  a  time  step  At  ths  simulated  molecules  move 
to  a  new  position,  and  the  macroscopic  velocity  and  temperature  of  a  cell  are  obtained  as  the  weighted 
averages  of  the  IP  velocity  m;  ^  and  IP  energy  e[  of  simulated  molecules: 


(4) 

k=\  k=\ 

^sm 

r=e'/(0.5<:)  =  (£e;/W.J/(0.5*:). 

(5) 

k=\ 


Here  is  the  total  number  of  the  simulated  molecules  in  the  cell,  T  and  e'  are  the  translational 
thermal  temperature  and  the  translational  kinetic  energy  based  on  the  ith  component  of  thermal 
velocity.  The  increment  of  the  IP  density  of  a  cell  during  At  is  (according  to  the  mass  conservation 
equation  (3)): 

An  =  -(At  /  AF) Jj  nU  ■  MS .  (6) 

The  macroscopic  density  is  updated  accordingly.  Meanwhile  the  increments  of  IP  values  u  and 
e‘  during  time  step  At  are  (according  to  the  momentum  and  energy  conservation  equations  (1),  (2)): 

Au  —  —(At  /(m«AF))JJ  pldS  (7) 

Ae’ = -(At /(mnAV)) q  -1(13,1=1,2,3  (8) 

where  AV  is  the  cell  volume.  The  surface  integrals  are  calculated  on  the  upper  and  lower  sides  of  each 
cell.  The  under  integral  values  in  (6),  (7),  (8)  are  the  number  flux,  normal  momentum  flux  and  heat 
flux,  and  are  calculated  from  the  assumption  that  the  molecules  in  each  cell  are  in  equilibrium  having 
the  cell’s  velocity,  density  and  temperature  as  their  macroscopic  characteristics.  For  example,  the  heat 
flux  across  the  boundary  of  cell  y=const  for  the  kinetic  energy  of  the  y-component  ( i  =2)  of  the  thermal 
velocity  is  (s  =  U l{2kTI ) 

q'f  =  +  /)exp(-5^)  ±  4^s(\.5  +  s^)(l  ±  erf(s))]  (9) 

whereas  that  for  the  energy  of  the  x  and  z  components  ( /  =  1,/  =  3 )  is 
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In  case  of  s=0  =kTr„,q‘"^  =q'y^  =0.5kTr„,  where r„  is  the  number  flux  in  a  stationary 

gas:  =  n(kT/(27im)y^^^^ .  The  difference  between  the  expressions  for  ciy^  and  {=q'y^)  is  the  reason 

why  the  kinetic  energies  for  different  components  of  thermal  velocity  are  introduced  separately.  When 
taking  the  surface  integral  on  each  side  of  the  cell,  the  energy  flux  consists  of  the  outlet  energy  flux 
from  the  cell  and  the  inlet  energy  flux  from  the  neighboring  cell,  and  the  pressure  p  is  obtained  as  the 
average  of  the  pressure  of  the  cell  and  the  pressure  of  neighboring  cell.  The  increments  Au  and  Ae'  are 
added  to  the  IP  values  of  the  simulated  molecules  in  the  cell,  the  updated  values  are  used  for  the  next 
calculation. 

During  collision  common  post-collision  IP  velocities  are  calculated  as  the  mass  weighted  averages 
of  the  pre-collision  values  and  are  assigned  to  each  of  the  encountering  molecules,  and  equipartition  of 
energy  is  applied  to  the  distribution  of  the  post  collision  IP  energies.  When  reflecting  from  walls,  for 
specular  reflection  IP  velocity  changes  its  sign  of  the  normal  to  the  wall  component  only,  and  for 
diffuse  reflection  it  gains  the  wall’s  velocity;  and  for  a  part  of  all  molecules  the  IP  energy  gains  the 
values  corresponding  to  the  wall  temperature  and  for  the  rest  (1-a)  part  the  IP  energy  remains 
unchanged,  where  a  is  the  accommodation  coefficient.  The  procedure  is  iterated  xmtil  convergence  is 
reached.  On  the  surface,  the  energy  exchange  is  implemented  by  concrete  counting  the  energy  carried 
by  the  simulated  molecules,  so  on  the  boundary  the  energy  flux  on  the  boundary  side  of  the  integral  in 
equation  (8)  is  not  calculated  repeatedly. 

As  we  have  seen,  in  the  equilibrium  state,  the  heat  flux  across  a  surface  with  normal  of  i  =  2  is 
((0.5A:T"‘  +kT"^  +0.5kT"^)r^,  so  a  molecule  when  going  through  such  a  surface  carries  in  average  an 
energy  of  0.5k(T‘"' +  2T‘"^) ,  and  the  heat  exchange  of  molecule  with  the  wall  surface  is 

=  Q.5k{ATj"^  +  AtJ"^  +2AtJ"^).  During  the  sampling  time  t  ,  the  net  energy  exchange  E  is  the  sum  of 
the  heat  exchange  E^  of  all  molecules  participating  in  sampling.  The  heat  transfer  to  the  wall  is 
q=  EI(t-AA)  ,  where  AA  is  the  area  of  the  wall  surface. 

RESULTS  AND  DISCUSSION 

The  problem  of  rarified  gas  confined  between  two  parallel  plates  with  different  temperatures  is 
simulated  by  the  IP  technique  described  above  for  a  temperature  difference  of  AT ITq=  0.14  and  with 
accommodation  coefficients  of  the  upper  and  lower  plate  =  a_{=  a)  =  0.826  . 

The  table  I  shows  the  variation  of  normalized  heat  flux  with  Kn  in  the  problem,  where  is 

the  analytical  result  of  gjfor  A^«  =  oo  given  by  [5]  {Qp^^  =  -{2{k  I  ^  ^  pqa{2  -ay^  AT  I T^)* .  The 

present  result  is  compared  with  the  solution  of  the  Boltzmann  equation  the  result  of  the  moment 
method  and  the  calculation  by  DSMC  method  carried  out  in  this  paper.  When  the  Kn  number  is 
small,  the  heat  flux  of  IP  method  is  the  smallest  among  all  methods,  and  the  accuracy  is  not 
satisfactory.  However,  when  the  Kn  number  is  high,  the  heat  flux  of  IP  method  is  in  good  agreement 
with  DSMC  procedure,  and  we  can  anticipate  when  Knudsen  number  tends  to  infinity  (i.e.  the  flow  can 
be  regarded  as  collisionless  flow),  the  results  of  both  methods  will  tend  to  the  analytical  result  Qp^ . 
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Table  I.  The  normalized  heat  flow  Q,  !  Qfm  various  Knudsen  numbers.  (  AT  IT(^=  0.14  ) 


Kn 

Solution  of  Boltzmann 
equation^^^ 

0.0658 

0.2132 

0.3018 

0.2538 

0.2716 

0.1942 

0.4230 

0.5336 

0.4843 

0.5239 

0.7582 

0.7072 

0.7813 

0.7558 

0.8112 

10. 

0.9584 

0.9684 

100. 

0.9865 

0.9906 

The  density  distribution  in  the  gas  for  different  Kn  obtained  by  the  IP  procedure  is  shown  in  Fig.  1 
together  with  the  experimental  data  the  solution  of  the  non-linear  Boltzmann  equation^^^  and  the 
result  of  DSMC  simulation  obtained  in  this  paper.  Density  and  Y  are  normalized  by  the  centerline  value 
and  D,  respectively.  In  general,  the  calculation  of  IP  method  is  in  fair  agreement  with  the  results  of 
other  methods.  At  about  Kn=0.5,  the  departure  of  IP  result  from  DSMC  calculation  and  the  solution  of 
the  non-linear  Boltzmann  equation  is  the  largest. 

Figure  2.  shows  the  temperature  distribution  of  IP  method  with  the  solution  of  the  non-linear 
Boltzmann  equation  and  the  result  of  DSMC  simulation  (the  temperature  is  normalized  by  Tq).  The 
agreement  among  these  methods  is  fair. 

Figure  3.  shows  the  density  and  temperature  distribution  of  IP  technique  and  DSMC  method  for 
Kn=0.7582  with  =a_(=a)  =  0.826  and  dT/ Tq  =  0.03  .  When  the  temperature  differences  {AT !  T^) 
decrease,  the  agreement  of  density  and  temperature  distribution  between  IP  result  and  DSMC 
calculation  is  better. 

The  above  comparison  shows  that  the  present  IP  procedure  gives  satisfactory  results  for  both  heat 
transfer  and  density  and  temperature  distribution  for  Kn  numbers  near  the  free  molecule  limit.  As  for 
small  Kn  numbers,  it  can  yield  density  and  temperature  distribution  results  in  fairly  good  agreement 
with  DSMC  simulations,  but  can  not  provide  satisfactory  heat  transfer  results. 

Of  course,  the  IP  procedure  has  the  advantage  of  prompt  convergence  which  is  not  so  important  for 
the  present  case,  but  is  crucial  for  complicated  MEMS  flows.  In  this  case,  when  the  IP  sample  size  is 
about  4,000,000  per  cell,  the  convergent  result  (when  the  y-component  (i=2)  of  the  macroscopic 
velocity  is  less  than  0.002)  is  achieved;  meanwhile  DSMC  method  should  have  the  sample  size  of 
about  1,000,000,000  per  cell  to  reach  convergence. 

Unlike  the  isothermal  flow  cases  where  the  IP  simulation  results  yield  good  agreement  with  DSMC 
calculation  and  other  exact  solutions,  the  present  IP  technique  for  non-isothermal  flow  provides  results 
only  in  fairly  satisfactory  agreement  with  DSMC  simulation  and  other  solutions  and  for  the  heat  flux  of 
small  Kn  number,  the  agreement  is  not  satisfactory.  The  technique  deserves  further  development. 
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density  |  |  density  |  |  density 


(a) 


FIG.  3.  Density  and  temperature  distribution  for  Kn=0.7582  with  ad^O.826  and  AT  /  TJj  =  0.03  . 
Here  ,  - :  IP  result;  - :  DSMC  calculation. 
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